plot(c(-1,6),c(-1,10),axes=F,type='n',ylab='Kullback-Leibler scale',xlab='')
 box()
 axis(2,at=seq(0,10,2),labels=c('','',0,rep('',3)))
 abline(h=4,lty=5,lwd=2,col='seagreen')
 
 #truth
 segments(0,1,5,1,lwd=2,col=2)
 text(2.5,1,expression(italic(f) == 'Truth'),col=2,pos=1)

 #model g 
 segments(0.25,9,2.25,9,lwd=2,col=4)
 text(1.25,9,expression(paste('Model ',italic(g))),col=4, pos=3)
 arrows(2.1,9,2.1,1,code=3,length=.075,col=4,angle=45,lty=2)
 text(2.1,5,expression(italic(I)(italic(f),italic(g))),pos=2,col=4) 
 arrows(0.35,9,0.35,4,code=3,length=.075,col=4,angle=45)
 text(0.35,6.5,expression(AIC(italic(g))),pos=2,col=4)
 
 #model h
 segments(2.75,6,4.75,6,lwd=2,col='grey20')
 text(3.75,6,expression(paste('Model ',italic(h))), col='grey20',pos=3)
 arrows(2.9,6,2.9,1,code=3,length=.075,col='grey20', angle=45,lty=2)
 text(2.9,3.4,expression(italic(I)(italic(f), italic(h))), pos=4, col='grey20')
 arrows(4.65,6,4.65,4,code=3,length=.075,col='grey20',angle=45)
 text(4.65,5,expression(AIC(italic(h))),pos=4,col='grey20')